Large-N expansion based on the Hubbard-operator path integral representation and 

its application to the t — J model 
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In the present work we have developed a large-N expansion for the t — J model based on the path 
integral formulation for Hubbard-operators. Our large-N expansion formulation contains diagram- 
matic rules, in which the propagators and vertex are written in term of Hubbard operators. Using 
our large-N formulation we have calculated, for J = 0, the renormalized 0(1/N) boson propagator. 
We also have calculated the spin-spin and charge-charge correlation functions to leading order 1/N. 
We have compared our diagram technique and results with the existing ones in the literature. 
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I. INTRODUCTION 

The role of electronic correlations is an important and 
open problem in solid state physics. Its close connection 
with the phenomena of high-T c superconductivity!!! makes 
this problem relevant in present days. 

One of the most popular models in the context of high- 
Tc superconductivity is the t — J model. A natural rep- 
resentatioit,of the t — J model is in terms of Hubbard 
operatorsoH. Although the t — J model is quadratic, 
when is written using Hubbard-operatorsa, there are sev- 
eral difficulties in the calculation of physical quantities. 
These difficulties are mainly: a) the complicate commuta- 
tion rules of the Hubbard operatorscl (X-operators) and, 
b) any of the model parameters can be taken to be small 
in a perturbative evaluation of the observables. 

To handle the point a), a popular method is to use 
slave particles (slave boson and fermion) to decouple the 
original A-operator in usual bosons and fcrmionsa For 
point b), one of the nonperturbative techniques (which 
will be relevant in the present paper) is the large-A ex- 
pansion. N is the number of electronic degrees of freedom 
per site (see below) and 1 /N can be considered as a small 
parameter. 

The large-A expansion was developed in the frame- 
work of skive boson representationHfl and, using Baym- 
KadanoffLLp. functional theory in terms of the X- 
operatorsolj. 

In this paper we develop a large-TV expansion for the 
t — J model based on the path integral representation for 
A-operators. nn 

In a series of papersE30 we have studied the path 
integral formulation for Hubbard-operators. Our point 
of view is to write a path integral for t-J model by 
using directly the Hubbard A-operators as fields vari- 
ables, without any decoupling. Our starting point was 
the construction of a particular family of first-order con- 
strained LagraxLgians by using the Faddeev-Jackiw and 
Dirac mcthodsEj. We showed that the constrains of the 
model are second class. The canonical quantization of 



this constrained theory leads to the commutation rules 
of the Hubbard-operators. Next, by using path-integral 
techniques, the correlation functional and effective La- 
grangian were constructed. 

In Rcf . [ |ll[] , we found a particular family of constrained 
Lagrangians and showed that the corresponding path- 
integral can be mapped to that—of, the slave-boson rep- 
resentation in the radial gaugeE3tJ. The path-integral 
formulation defined in [ |ll[] is the starting point for our 
large-N expansion. 

In section II we give a summary of the formalism 
which will be useful for understanding our large-N ex- 
pansion. In section III, we develop our large-N method. 
We present explicit calculations for the case J = of 
the t — J model ([/-infinite Hubbard model). In section 
IV, we calculate spin-spin and charge-charge correlation 
functions to leading order of large-N expansion. Finally, 
in section V we give our main discussions and conclu- 
sions. 



II. PRELIMINARIES AND DEFINITIONS 

In this section we give a summary of the main re- 
sults obtained in [11], which are the starting point for 
our large-N expansion. 

As we mentioned in the introduction, although the t— J 
model is quadratic when is written using X-operators, 
the complicate commutation rules between them make 
this model highly nontrivial. To attack the problem, our 
method works as follows: 

a) We construct a Lagrangian written in terms of Hub- 
bard operators 



£ = 5> M/3 (A)A i "' 3 -H(X) 
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where the t — J model Hamiltonian H(X), with the 
addition of a chemical potential /i, reads as: 
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In (1) the coefficients cii Q ^(X) are unknown and must 
be determined. 

The indices a and /3 can take both value (empty 
state) or spin index a = ± (up and down state, respec- 
tively). The five Hubbard X-operators X aa and X 00 are 
boson-like and the faur Hubbard X-operators X a0 and 
X 0a are fermion-likea. In (2), Uj and Jij are hoping and 
exchange parameters, respectively, between sites i and j. 

b) We introduce (via Lagrangian multipliers A Q ^) a set 
of bosonic constrains VL a ^ which must be also determined. 

c) After these definitions and proposals, we impose 
that the Dirac bracketstj between X-variables obtained 
using our constrained theory verify the correct commuta- 
tion rules for Hubbard operators. With this prescription, 
we obtain a set of differential equations for the coefficients 
a>iap{X) and constrains Sl a . A particular solution for the 
coefficients is 
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The boson-like Lagrangian coefficients are all zero. 
The constrains are second class because none of them 
commutes with all the other ones. 

Consequently, the dynamics is given by the Lagrangian 



L(X,X) = - l -J2 { -^^ 
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and the constrains (4) and (5). 

Now, in order to write the path integral representa- 
tion we use the method of Faddeev-SenjanovichEj. This 
method is currently used in quantum field theory to ob- 
tain the path integral of constrained systems. 

The partition function reads as: 



z = [ vxf s[x°° + x r - 1] t[ x r' - 

J a 

x (sdetM AB )Jexp (i J dtL(X,X)) . 
where the superdeterminant 
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[sdetMABY = 1/ 
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is just the superdeterminant of the Dirac matrixO (or 
equivalently the sympletic matrix of the Faddeev-Jackiw 
method) formed with the set of all the second class con- 
strains of the theory. 

The power two on X 00 in (sdetMAB)^ appears because 
we are working with two spin projections. This point will 
be very important for our large-N extension in the next 
section. 

Note the first constrain (4) (the first <5-function in (7)) 
is the completeness condition which means that "double 
occupancy" at each site is forbidden. The constrains (5) 
(the second <5-function in (7)) is obtained in our formal- 
ism as a consequence of imposing that the commutation 
rules of X-operators must be fulfilled. 

From the above equations our theory corresponds to a 
configuration in which the bosons are totally constrained 
and the dynamics is carried out only by the fermions. 
As it was shown in Ref[ll], the path integral formalism 
corresponding to this dynamical situation is mapped in 
the slave boson representation. 

Equation (7) is our main formula to develop the large- 
N expansion. 



III. THE LARGE- TV EXPANSION 

In this section, we present the large- TV expansion in the 
framework of the path integral representation for Hub- 
bard operators. We present our explicit results for the 
J = case (L^-infmite Hubbard model) which will be very 
useful to show how our method works and to compare our 
results with previously ones found in the literature. 

Our starting point is the partition function (7) written 
in the Euclidean form: 

X a ® X® a ' 



z = J vxf S[X°° + ]T Xf ° l] 5[xf Y1 ,„ 

H{X) (6) x (sdetM AB )f exp {- [ dr L E (X,X)) . (9) 



The Euclidean Lagrangian Le[X,X) in (9) is: 



First, we make the following steps in the path integral 
(9): 

a) we integrate over the boson variables X 00 using the 
second (5-function in (9), 

b) the spin index a = ±, is extended to a new index p 
running from 1 to N. In order to get a finite theory in 
the X-infinite limit, we re-scale the hoping ty to tij/N 
in the Hamiltonian. 
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c) the completeness condition (Xf° + J2 P x f P = N l^) 
can be exponentiated, as usual, by using the Lagrangian 
multipliers Aj, 

d) in order to prepare the path integral for the large-N 
expansion, we write the boson fields in terms of static 
mean-field values and dynamic fluctuations 



Xf i = Nr {\ + SR i ) 



(11) 



A, = A + SX t , (12) 
e) finally, we make the following change of variables 



f + = — 



fil 



xf, 



(13) 



(14) 



By following the steps a-e, we find the effective La- 
grangian for the J = case 

Leff = ~2 E {jipfip + fvpfip) (1 _|_ SR.) 

i,p ^ %l 

N 1 
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In (15), L ghost is defined via the exponentiation of the 
superdenterminant using ghost fields. After the exten- 
sion to large-TV, the superdeterminant becomes: 



(sdetM AB )f = 



(Nr ) 



N 



(16) 



The numerator of (16) together with the Jacobian of the 
transformations (11-14) contribute to the path-integral 
normalization factor. The denominator of (16) can be 
seen as a N x N diagonal matrix and the integral repre- 
sentation is given in terms of complex boson ghost fields 
Z p . Therefore, 

L ghost (z) = - E 4 (tTsr) Zip ■ (17) 

ip ^ L ' 

To implement the 1/N expansion, the nonpolynomial 
L e ff should be developed in powers of SR. Up to order 
1/N, it is sufficient to retain terms up to SR 2 . Then, 
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(18) 
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where we have changed fi to \l — A and dropped con- 
stant and linear terms in the fields. 

Looking at the effective Lagrangian (18), the Feynman 
rules can be obtained as usual. The bilinear parts give 
rise to the propagators and the remaining pieces are rep- 
resented by vertices. Besides, we assume the equation 
(18) written in the momentum space once the Fourier 
transformation was performed. 

In summary, the Feynman rules are: 

i) Propagators: We associate with the two component 
SX a = (SR , SX) boson field, the propagator 



D(0)ab{<l,U n ) - 
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which is represented by a dashed line in Fig.l connecting 
two generic components a and b. 

The quantities q and u) n are the momentum and the 
Bose Matsubara frequency of the bosonic field, respec- 
tively. 

We associate with the N-component fermion field f p , 
the propagator 



G 



(O)pp' 



(k,v n ) = - 



iv n - (e k - fi) 



(20) 



which is represented by a solid line in Fig.l connecting 
two generic components p and p'. 

In (20) the electron bare dispersion was defined as Ek — 
—2tr (cosk x + cosky), where t is the hopping between 
nearest neighbors sites on the square lattice. 

The quantities k and v n are the momentum and the 
Fcrmionic Matsubara frequency of the fermionic field, re- 
spectively. 

We associate with the N-components ghost field Z p i, 
the propagator 



Sppf 



(21) 



which is represented by a dotted line in Fig.l connecting 
two generic components p and p\ 

ii)Vertices: The expressions of the different three-leg 
and four-leg vertices are 



A^ = (-l) Q(„ + ,/)+„; l) 



(22) 
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representing the interaction between two fermions and 
one boson (Fig. la); 

A# = (-1) ( " 7 ) K (23) 

representing the interaction between two fermions and 
two bosons (Fig. lb); 



Kp> = (-!)(V . 0), 



(24) 



representing the interaction between two ghosts and one 
boson (Fig.lc); and 



-1) 



-1 




3pp' 



(25) 



representing the interaction between two bosons and two 
ghosts (Fig. Id). Each vertex conserves the momentum 
and energy. 

In the iV-infinitc limit (order O(l)), we only have free 
fermions with renormalized band due to correlations, 
e k = —2tro(cosk x + cosky). For a given value of fi, r$ 
(which from (11) and the completeness condition is equal 
to 6/2, where (5 is the hole doping away from half- filling) 
must be determined self-consistent]*- n This result is in 
agreement with previous calculatiorJI3'El. 

The bare boson propagator (19), which is of order 
0(1/N), is renormalized by a series of electronic bubbles 
which contribute also in 0(1/N). Then, in 0(1/N), the 
boson propagator is dressed by electronic interactions. 

By looking at the diagrammatic, it can be seen that the 
irreducible boson self-energy Tl a i, is given by the sum of 
the contributions corresponding to the one-loop diagrams 
(in order 1/N) seen in Fig.le. 




(c) 



(d) 
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...Q... 



FIG. 1. Types of vertex, up to order 0(1/N), between two 
fermions and one boson (a), two fermions and two bosons 
(b), two ghosts and one boson (c) and, two bosons and two 
ghosts (d). (e) represents the sum of all one- loop diagrams 
contributing to the irreducible boson-self energy. 
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[n F {e k+q - fi) - n F (e k - n)] 
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and, 



N 



1 V V — ^ 

Uxx(q,ui n ) = -— 



[n F (e k+q - n)- n F (e k - fi)] 
—iuj n + e k+q — e k 



(28) 



From the Dyson equation (D ab ) 1 = (-D(o)ab) 1 - U ab 
the dressed components D ab of the boson propagator can 
be found. 

Our expressions (26)-(28) agree with the previously ir- 
reducible boson self energies calculated in the context of 
the large-iV" slave boson approachd. 

The ghost fields interact only with the boson fields 
(Fig.lc and Fig. Id). The only role of the ghost fields, in 
order 1/N, is to cancel infinities in the boson self energy 
II a b due to the frequency dependence of our vertices (22) 
and (23). For the evaluation of quantities to higher orders 
than 1 /N, we must include and check the contribution of 
the ghost fields in each particular calculation. Moreover 
it should be noted that in this case we must consider 
terms beyond SB 2 in developing the nonpolinomial parts 
of L e ff. Although this seems to be complicated, all the 
contributions can be controlled by the small number 1 /N 
via the vertices and propagators. At this point we want 
to remark that so far not many calculations go beyond 
the order 1/N. In order 1/N, diagramatic rules are well 
defined in the present paper. 

In summary, we have developed a diagrammatic tech- 
nique for a large- N expansion in the same spirit of the 
large- N expansion in quantum field theory It means, via 
the order of the propagators and vertices, we can deter- 
mine the order of the diagram contribution. Note that 
our Green's functions are calculated in terms of the origi- 
nal Hubbard operators. For example, from (13) and (14) 
we can see that our fermions fi p are always proportional 
to the Fermi-like X-operator X op and not only— in the 
leading order like in the slave-boson formulationt^l. 



The presence of the last two diagrams in Fig.le involv- 
ing ghost fields is very important. It is possible to show 
that these two diagrams exactly cancel the infinities, due 
to the frequency dependence of our vertices, of the two 
first diagrams in Fig.le. Therefore, using our Feynman 
rules, the final (and finite) result for each component of 
the boson-self energy H a b is: 



IV. CHARGE-CHARGE AND SPIN-SPIN 
CORRELATION FUNCTION 

In this section we calculate numerically the charge- 
charge and spin-spin correlation functions on the square 
lattice for the nearest neighbor hopping t. We choose 
t = 1.0 as the energy unit and the temperature T = 0.0A". 
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We can define retarded .density-density D and spin- 
spin S Green's functions asEJ 



D ij = ^Y,< T r x r x r> 



(29) 



and 



~S l3 =< T T Xf"Xf >, 



(30) 



respectively. 



Using V Xf q = N/2 - X°° and (11) we find for D, 



0(1), in the Fourier space 



D(Q,u n ) = -N(-) Drr(Q, uj n ) 



(31) 



Using the constrains (5), it is easy to prove that S, in 
0(1), is the electron bubble formed with two bare elec- 
tron Green's functions (20). The analytical expression 
for S agrees with U\\. 

In Fig. 2, we show D(Q,u>) = — Im[D(Q, u>)] and 
S(Q,uj) = — Im[S(Q, (J)] using rj = O.lt in the analyt- 
ical continuation iu> n = u> + ir). 

To make explicit comparison with 
exact diagonalizationE£l and analyticaO calculation wc 
choose the same densities and momenta Q of these two 
references. 

In Fig. 2a, we plot D(Q,oj) (solid line) and S(Q,u>) 
(dashed line) for doping 6 — 0.12. In Fig. 2b, we present 
similar results for doping 8 — 0.78. In agreement 
with analytical calculation based on A-operatcp,canoni- 
cal approached (see also slave boson calculationEj), there 
is no collective excitations (like magnons) in the spin-spin 
correlation function to leading order of the 1/N expan- 
sion. The spin-spin correlation functions is determined 
by the particle-hole bubble with the renormalized elec- 
tron dispersion by correlations. In contrast to S(Q 1 uj) 1 
D(Q,lo) contains collective excitations (zetp sound) that 
can live inside or outside of the spectrumll3. 

For low doping (Fig. 2a) the collective excitations lie 
outside the particle-hole spectrum. Therefore, they show 
up as sharp peaks. In contrast, for high doping (Fig. 
2b) the collective excitations lie inside the particle-hole 
spectrum and they show up as broad peaks. The domi- 
nance of collective excitations in D, for low doping, can 
be clearly seen for Q — (tt,tt). 

It is also very important to point out that spin and 
charge fluctuations occur at low doping (Fig. 2a), near 
half filling, on different energy scales. This picture dis- 
appears at high doping as it can be seen in Fig. 2b. This 
behavior is not present in the uncorrelated case where 
spin and charge excitations are essentially degenerate. 

The agreement between our results and previous ones 
(see Fig.l of Ref[Q) is excelent, proving the consistency 
of our large- TV method. 
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FIG. 2. Spin-spin (dashed lines) and charge-charge (solid 
lines) correlation functions as a function of uj for several mo- 
menta for dopings 8 = 0.12 (a) and <5 = 0.78 (b). 



V. CONCLUSIONS AND DISCUSSIONS 

Based on our Hubbard operator path-integral repre- 
sentation, defined for the first time in Ref. [ |ll| , we have 
developed a large- N expansion for the t — J model. 

Because our path integral is written in terms of X- 
operators, we do not need any decoupling scheme for the 
Hubbard operators. Then, we do not invoke a gauge 
invariance ot,Bose condensation, like in the slave boson 
formulationEj. Our dynamical variables are the original 
Hubbard operators and, there is no need to introduce 
spinons and holons excitations. 

The use of a path integral allows us to define a Feyn- 
man diagram technique which is a very importapt and 
interesting problem in electron correlated systemsO. 

As we mentioned previously, the large-A expansion 
was formulated using canonical quantization in terms of 
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X-operatorB'EI or slave boson path integral approach!^. 
The first one cannot be defined in a diagrammatic form 
while the second one can be developed in terms of Feyn- 
man diagrams but needs the introduction of spinons and 
holons. In the formulation presented in this paper we 
have diagrammatic rules: propagators and vertices in 
terms of the original X-operators. 

One important point is that in our diagram technique 
the fermion operator f a is proportional to the fermion 
X- operator X 0(T for all orders in the large N expansion. 
This is very important for the interpretation of the ex- 
perimental spectral susceptibilities. In the slave boson 
approach the spinon is only proportional to the fcrmipn 
X 0a operator in the leading order of 1/N expansionll3. 

It is also important to note that our diagram tech- 
nique contains ghost fields that cancel infinities to order 
0(1/ N) in the boson propagator. In the slave boson ap- 
proach in the radial gauge, the ghost fields necessary to 
treat the Jacobian, cancel infinities in 0(1/N 2 )t3. 

Using our large- N expansion, we have calculated the 
boson propagator and the charge-charge and spin-spin 
correlation functions in 0(1). In sections III and IV 
we have shown that these— quantities agree with those 
calculated-iby slave bosonE^I and A-operator canonical 
approachHl. These results, which are far nontrivial, show 
the consistency of our formulation. 

We showed that in the N-infinitc limit we have well de- 
fined electronic excitations with renormalized band due 
to correlations. This means (and we think) that the large- 
N expansion, in the form presented in this paper, is a 
better approximation to describe high-Tc cuprates in op- 
timal doping and overdoped than in underdoped close to 
antiferromagnetism. 

We believe that our formulation will be useful (and 
simpler than previous ones) to calculate quantities in 
the next order in 0(1/N), such as superconductivity and 
electronic self-energy corrections. We are also planning 
to extend the method to the case of finite J. 
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